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ABSTRACT 

Aims. We intended to measure the radial velocity curve of the supergiant companion to the eclipsing high mass X-ray 
binary pulsar EXO 1722-363 and hence determine the stellar masses of the components. 

Methods. We used a set of archival K s -band infrared spectra of the counterpart to EXO 1722-363 obtained using ISAAC 
on the VLT, and cross-correlated them in order to measure the radial velocity of the star. 

Results. The resulting radial velocity curve has a semi-amplitude of 24.5 ± 5.0 km s _1 . When combined with other 
measured parameters of the system, this yields masses in the range 1.5 ± 0.4 - 1.6 ± 0.4 M for the neutron star and 
13.6 ± 1.6 - 15.2 ± 1.9 Mg for the B0-1 la supergiant companion. These lower and upper limits were obtained under 
the assumption that the system is viewed edge-on (i = 90°) for the lower limit and the supergiant fills its Roche lobe 
(/3 = 1) for the upper limit respectively. The system inclination is constrained to i > 75° and the Roche lobe-filling 
factor of the supergiant is /3 > 0.9. Additionally we were able to further constrain our distance determination to be 7.1 
< d < 7.9 kpc for EXO 1722-363. The X-ray luminosity for this distance range is 4.7 x 10 35 < L x < 9.2 x 10 36 erg 
s _1 . 

Conclusions. EXO 1722-363 therefore becomes the seventh of the ten known eclipsing X-ray binary pulsars for which 
a dynamical neutron star mass solution has been determined. Additionally EXO 1722-363 is the first such system to 
have a neutron star mass measurement made utilising near-infrared spectroscopy. 

Key words. binaries:eclipsing - binaries:general - X-rays:binaries - stars:individual:EXO 1722-363 

1. Introduction Observations made of EXO 1722-363 (alternatively des- 

, . ignated ICR J17252-3616), in 1987 by the Ginga X-ray 

The precise form of the neutron star (NS) equation of state satemte were the first to detect pu l sat ions. These pulsa- 

is still unknown. Despite much theoretical work aimed at tiong were found to haye a 4 i3. 9± o.2 s period (jTawara et all 

determining this fundamental aspect of astrophysics, to rj 989 ). ver an 8 hour period the source appeared to 

eliminate some of the contending theories we must turn to yary substantially in fluX; decreasing from 2 mcrab to 

observational data. Presently the only means of determm- q.2 - 0.3 mcrab within the 6-21 keV band. It was sub- 

mg the mass of neutron stars m accretion driven systems is seque ntly found that the X-ray flux remained persistent 

by observing eclipsing X-ray binary pulsars. Unfortunately from 2Q _ 6Q keVj however a cutoff point was found 

only 10 such systems are currently known, and only 6 of above this flux leyel in which the source was unde tectable. 

these have previously had mass measurements made (e.g. EX 1722-363 at maximum flux and assuming a distance 

Ash etaJJ 019991), i Quamtrell et al.| d2Q03j),| Val Baker et al.| f 1Q kpCi had a luminosity calculated as 5 x 10 36 erg s" 1 



(20050). In this paper we present the preliminary results (| T awara et al.lflllh . Later observations by the Rossi X- 



from our on-going work on the mass of the High Mass X- my Tlming Exp i orer ( RX TE) revealed the eclipsing na- 

ray Binary (HMXB) accretion driven pulsar EXO 1772- ture of thig syst e m with the eclipse duration determined 

363. The counterpart star within this HMXB is heavily ob- as L7± .! days (|Corbet et al.ll2Tjol . Subsequent observa- 

scured and reddened, necessitating for the first time the tiong by INTEGRAL followed up by XMM-Newton in 2004 

utilisation of near-infrared spectroscopy to construct the led to a further refinement of the spin and orbital periods 

Radial Velocity (RV) curve m order to obtain an accurate to 413.851 ± 0.004 s and 9.7403 ± 0.0004 days respectively 

mass solution. (|Thompson et al.ll200l . 

* Based on observations carried out at the European Southern XMM-Newton observations allowed the source posi- 

Observatory under programme ID 077.B-0872(A) tion to be determined more precisely (with an uncer- 
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tainty of 4") at RA(2000.0) = 17 /l 25 m 11.4 s and Dec = 
-36°16'58.6". EXO 1722-363 lies within the Galactic plane 
and as it is heavily reddened, unsurprisingly the counter- 
part star could not be detected optically. An infrared coun- 
terpart was found lying 1 " from the X-ray source position 



(jZurita Heras et al 



ying r 
, 2006) 



with a corresponding entry in 



the 2MASS catalogue, 2MASS J17251139-3616575 (with 
JHK magnitudes J = 14.2, H = 11.8 and K s = 10.7). 
Examination of near infrared K-band spectra obtained with 
the ESO ISAAC instrument led to our determination of 
the spectral classifi cation of the mass donor as BO - Bl la 
(jMason et al.ll2009l ). 



2. Observations and Data Reduction 

In our previous work, we only had single epoch K s -band 
spectra of the mass donor in EXO 1722-363, but in order 
to determine a dynamical mass solution, radial velocities at 
a range of orbital phases are required. Fortunately we were 
able to locate a series of K s -band spectra held within the 
ESO ArchiveQ which were obtained over 26 nights between 
24th May, 2006 and 4th August, 2006. 26 pairs of spec- 
tra were centred on 2.1 fim and 26 pairs centred on 2.2 fini. 
No observations of radial velocity standards appear to have 
been taken with any of the science spectra, however there 
are telluric standards available to enable removal of atmo- 
spheric features from the target spectra. 

The data were reduced using the ISAAC pipeline in 
conjunction with the data browsing tool GASGANO 0. 
Unfortunately within the 2.2 fim EXO 1722-363 dataset, 
the counterpart to EXO 1722-363 had been incorrectly 
identified and the telescope had been mis-pointed, so the 
spectra were unusable. Additionally only 11 of the 26 spec- 
tra from the dataset centred on 2.1 /xm turned out to be 
of sufficient quality to derive a radial velocity measurement 
(See Tabled]). The usable spectra were made in the SW 
MRes mode with a 0.6" slit. Resulting in spectra with a 
high S/N and resolution (R w 4200). The integration time 
for each pair of spectra was 700 s, with the resulting data 
having a count rate below 10 000 ADU; therefore no cor- 
rection for non-linearity was necessary. 

Flatfields were reduced using the pipeline recipe 
isaac_spc_flat and combined to produce a master flatfield. 
The wavelength calibration and ISAAC slit curvature dis- 
tortion was computed using OH skylines using the pipeline 
recipe isaac_spc_arc. Spectra produced by ISAAC have a 
high degree of curvature, to remove this the pipeline recipe 
isaac_spc^startrace computes the spectra curvature using 
both images and spectra of a star moving across the slit. 
Science spectra were obtained using the nodding technique; 
unfortunately only one nod was performed during the orig- 
inal observation, which we believe to be less than optimal. 
The two nodded science frames were then reduced using 
the products of the pipeline calibration recipes to produce 
a final reduced science spectrum. This process was then re- 
peated for each telluric standard. Telluric correction was 
then made using the standards shown in Table [TJ All spec- 
tra were reduced using standard IRAjQ routines; Fig [T] 




Wavelength (Microns) 

Fig. 1. Continuum normalised K s -band spectra centred on 
2.1 /itm of EXO 1722-363 in order of date from top to bot- 
tom. 



shows the stacked continuum normalised spectra ordered 
by date of observation. 



3. Data Analysis 

Radial velocities were determined by cross-correlating the 
region around the Hel 2.112/3 /im absorption line in each of 
the 11 archive spectra against the high signal-to-noise K s - 
band spectrum of EXO 1722-363, which we h ad previously 
obtai ned for spectral classification purposes (jMason et alj 
I2009T ) . The resulting velocities were then corrected to the 
solar system barycentre and are reported in Table [1] The 
spectrum highlighted in bold is that used as the reference 
spectrum (our previously obtained high signal-to-nois e K s - 
band spectrum of EXO 1722-363 (jMason et alj 12009)). In 
obtaining these final velocities, we determined the absolute 
velocity of the reference spectrum by fitting the positions 
of its absorption lines. 

From X-ray data there is no evidence that 
EXO 1722-363 has an ything other than a circular orbit 
(jThompson et al.1 [2007) , so we fitted the radial velocities of 
the supergia nt star with a simple si nusoidal solution. The 
ephemeris of Thompso n et al.l (2007) specifies the epoch of 
mid-eclipse as 

T(HJD) 



53761.68(4) + 9.7403(4)A>" 



(1) 



where N is the cycle number and uncertainties in brack- 
ets refer to the last decimal place quoted. At the epoch of 
our observations, the accumulated uncertainty in phase is 
formally only ~ 0.005, but nonetheless we fitted our data 
with two models, in one of which the zero phase was a free 
parameter and in the other of which it was not. 

Fitting our data with a sinusoid with just two free pa- 
rameters (RV amplitude and systemic velocity) yielded an 
amplitude of 17.6 ± 7.7 km s _1 and a systemic velocity 
of -6.5 ± 5.6 km s . In order to achieve a reduced chi- 
— squared of unity, the uncertainties on each RV data point 

1 I http : //archive . eso . o rg/eso/e so"_archive_main . html | had to be scaled to ± 17.5 km s" 1 . I n comparison, fitting 

2 http://ww.eso.org/sci/data-processing/software/pipel f ree parameters (i.e. 

3 http : //www . eso . org/ sci/data-processing/ software/ gasgano/. 

4 IRAF is distributed by the National Optical Astronomy Universities for Research in Astronomy, Inc., under cooperative 
Observatory, which is operated by the Association of agreement with the National Science Foundation. 
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Table 1. The phase, radial velocity and telluric standard for each EXO 1722-363 archive spectrum. 



Mid-point of Observations (UT) 


HJD 


Phase 


Radial velocity / km s 1 


Telluric Std 


Telluric Spec. Type 


2006 May 24.340 


245 3879.83 


0.130 


-18.61 ± 11.8 


Hip 093225 


B4V 


2006 Jun 22.025 


245 3908.53 


0.077 


-11.74 ± 11.8 


Hip 088109 


B5II 


2006 Jun 29.233 


245 3915.74 


0.817 


-21.39 ± 11.8 


Hip 070148 


B8III 


2006 Jul 14.167 


245 3930.67 


0.350 


14.92 ± 11.8 


Hip 089960 


B6V 


2006 Jul 17.248 


245 3933.75 


0.666 


-29.22 ± 11.8 


Hip 087616 


B9IV/V 


2006 Jul 20.122 


245 3936.63 


0.961 


-14.95 ± 11.8 


Hip 090336 


B7III 


2006 Jul 22.112 


245 3938.62 


0.165 


-10.09 ± 11.8 


Hip 094859 


B5V 


2006 Jul 24.115 


245 3940.62 


0.371 


21.21 ± 11.8 


Hip 090336 


B7III 


2006 Jul 25.105 


245 3941.61 


0.473 


6.89 ± 11.8 


Hip 089960 


B6V 


2006 Aug 01.111 


245 3948.57 


0.188 


-3.83 ± 11.8 


Hip 085548 


B9II 


2006 Aug 02.113 


245 3949.62 


0.295 


28.06 ± 11.8 


Hip 085548 


B9II 


2008 May 17.115 


245 4603.62 


0.439 


30.95 ± 11.8 


Hip 085008 


B5V 



with the addition of zero phase as a free parameter) gave 
an amplitude of 24.6 ± 5.0 km s _1 , a systemic velocity of 
—6.5 ±3.8 km s _1 and a phase shift of -0.13 ± 0.03. In this 
case, a reduced chi-squared of unity was achieved by scaling 
the uncertainties on each point to ± 11.8 km s _1 . Although 
the best-fit phase offset is discrepant with the accumulated 
phase uncertainty of the ephemeris, we prefer this fit and 
use the data from it subsequently; both fits are shown in 
Figure [2] in which the value for Ko is that resulting from 
fitting the radial velocities including a phase shift. 



dition, for circular orbits, 




Fig. 2. Radial velocity data for the supergiant star in 
EXO 1722-363. The solid line is the best fitting si- 
nusoid with three free parameters, the dashed line is 
that with a fixed zero phase in line with the published 
ep hemeris. The orb i tal ph ase is based upon the ephemeris 
of iThompson et all (|2007h . 



The masses of the system components may be deter- 
mined as follows. The mass ratio of the system q is equal 
to the ratio of the semi-amplitudes of the radial velocities 
for each star 

q=^ = ^ (2) 
q Mo K x { ' 

where Mx and Mo are the masses of the neutron star and 
supergiant star respectively, and Kx and Ko are the corre- 
sponding semi-amplitudes of their radial velocities. In ad- 



and similarly 



Mo = Ax P , (1 + g) 2 
27rGsin 3 i H> 



„ X = J^L( 1 + I 

27rGsnTi V q 



(3) 



(4) 



where i is the inclination to the plane of the sky and P 
is the orbital period. For EXO 1722-363, X-ray pulse tim- 
i ng delays yield the va lue of Kx as 226.1 ± 6.7 km s _1 
([Thompson et al.1 [2007). A value for the system inclination 
can be found from the geometric relation 



1/2 



cos 9 r 



(5) 



where 6 C is the eclipse half-angle, i?L is the Roche lobe 
radius of the supergiant, (3 is the ratio the supergiant 's 
radius to that of its Roche lobe and a is the separation 
between the centres of mass of the two stars. The Roche 
lobe radius may be approximated by 



— » A + B\ogq + C\og 2 q 
a 



(0) 



where the constants have been determined by 
Uoss fc Rappaportl (|1984l ) as 



A k 0.398 - 0.026S1 2 + 0.004ft 3 
B « -0.264 + 0.052ft 2 - 0.015ft 3 
C« -0.023 - 0.005ft 2 



(7) 

(8) 
(9) 



ft is the ratio of the spin period of the supergiant to its 
orbital period. For EXO 1722-363 we have assumed that 
the supergiant is close to Roche lobe-filling and is rotating 
synchronously with the orbit, so ft = 1 (although we note 
this assumption may not be entirely correct in this case), 
and the eclipse half a ngle is measured us ing RXTE data to 
be O = 31.8° ± 1.8° dCorbet et al.lf2005h . 

Hence the above set of equations allow the masses of the 
two stars to be determined in two limits. First, assuming 
that the supergiant fills its Roche lobe (in which case f3 = 1) 
we can find a lower limit to the system inclination i and 
upper limits to the stellar masses. Secondly, assuming that 
the system is viewed edge-on (in which case i = 90°) we 
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Table 2. System parameters for EXO 1722-363. 



Parameter 


Value 


Ref. 


Observed 






ax sin i / It sec 


101 ± 1 


[1] 


P 1 d 


9.7403 ± 0.0004 


[1] 


T go / HJD 


53761.68 ± 0.04 


[1] 


e 


< 0.19 


r 1 1 

w 


C / deg 


31.8 ± 1.8 


[2] 


Kq 1 km s" 1 


24.5 ± 5.0 


[3] 


Assumed 




n 


= 1 




Inferred 






K x 1 km s" 1 


226.1 ±6.7 




1 


0.107 ±0.022 





P 


1.000 


0.916 ±0.047 


i 1 deg 


75.2 ± 4.6 


90.0 


M x / M 


1.63 ± 0.38 


1.46 ± 0.38 


M / M 


15.2 ± 1.9 


13.6 ± 1.6 


a 1 R 


49.1 ± 9.1 


47.3 ±8.8 


Rl 1 R© 


28.0 ± 5.3 


27.0 ± 5.0 


-Ro / R© 


28.0 ± 5.3 


24.7 ± 4.7 



[1] Thompson et al. 2007; [2] Corbet et al. 2005 
[3] this paper 



can find a lower limit to the Roche lobe filling factor (3 
and lower limits to the stellar masses. Unfortunately, the 
spectra we have obtained from the ESO archive are not of 
sufficient quality to conduct a non-LTE model atmosphere 
analysis and hence make an accurate determination of the 
stellar radius which would break this degeneracy. 

In order to propagate the uncertainties in each param- 
eter, we performed a Monte Carlo analysis of the above 
equations to determine the system masses. The results in 
each limit are shown in Table [2] , and of course masses ly- 
ing between the extremes are also valid and correspond to 
values of i and f3 between their extremes. 

4. Discussion 

Although the 11 spectra reported here are of relatively low 
quality, and few in number, they still allow us to make a 
preliminary determination of the orbit of the supergiant 
in EXO 1722-363 and make a first measurement of the 
dynamical masses of the stellar components. The results 
are encouraging for a number of reasons. First, the re- 
sulting neutron star mass is consistent with the canon- 
ical mass of 1.4 M Q measured in most other eclipsing 
HMXBs, except for that in Vela X-l, (|Quaintrell et al.l 
2003). Second, the measured mass and radius of the su- 
pergiant, M ~ 13 - 15 Mq and R ~ 25 - 28 R Q , support 
the B0-1 la spectral classificati on that we have previously 
determined iMason et al.l 12009). This is illustrated by the 
Hcrtzsprung- Russell diagram plotted in Fig. which shows 
a close correspondence between the syst em primary and th e 
properties of other galactic field BSGs (jSearle et al.l l2008). 
While the similarity in temperature is to be expected - the 
value for the primary was adopted on the basis of its spec- 
tr al type, which in tur n has been calibrated by the analysis 
of ISearle et all (|2008l ). the radii for the field BSGs were 
determined via non-LTE model atmosphere analysis, while 
that for the primary is instead determined dynamically. 

The measurement of the stellar radius and hence bolo- 
metric luminosity, has in turn has allowed a more precise 
determination of the distance to the system by comparison 



to its observed photometric magnitude and reddening. The 
refined distance to EXO 1722-363 of 7.1 - 7.9 kpc results 
in an X-ray luminosity ranging from Lx min = 0.47 x 10 36 
to Lx max = 9.2 x 10 36 erg s -1 . However, due to the nature 
of the archive observations used for this work, the uncer- 
tainties on the mass and radius parameters are still rather 
large; it is our intention in the near future to propose and 
obtain more accurate VLT/ISAAC observations to further 
constrain the orbital solution parameters for this HMXB 
system. 

Finally, comparison to evolutionary tracks in Fig. [3] 
might suggest the primary had an initial progenitor mass 
of ~ 35 - 40Mq and hence the neutron star originated 
in a more massive star. However we caution that the bi- 
nary is highly likely to have undergone at least one episode 
of mass transfer in the past, rendering such a conclusion 
highly uncertain. As an exemplar we cite GX 301-2, a 
HMXB composed of a NS an d a B hypergiant w ith a spec- 
troscopic mass of 43±10M (n (|Kaper et alj 120061 ). However 
(jWellstein k Langerl fl999) propose a formation scenario in 
which two stars of comparable initial masses evolved via 
quasi conservative mass transfer into the current configura- 
tion post supernova; hence determining progenitor masses 
for both primary and neutron star based on the current 
system parameters is non trivial. 
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- Searle Galactic BSGs 

- Trundle SMC BSGs 

= Crowther Galactic BSGs 

- Repolust Galactic stars 
= EXO 1722 - 363 




Fig. 3. Position of EXO 1722-363 on the Hertzsprung- 
Russell (jSearle et al.l [2008) diagram alongside a sample 
of O and B sup ergiants from di ff ering locations, Galactic 



B supe rgiants, 



(2006 :), SMC B supe r- 



. ICrowther et alj - _ 
giants (iTrundle et alJ (120041): iTrundle fc Len"non1 (2005)) 
and Galactic O stars, iRepolust et al.1 (|2004l ). These are 



overplotted together with solar met allicity evolutionary 
tracks from iMevnet fc Maederl (|2000T ). Also shown is the 
lower and upper limits on the luminosity of EXO 1772-363. 



